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I. INTRODUCTION 

Landau's Fermi liquid theory^ is an immensely suc- 
cessful and widely used framework for understanding 
the properties of interacting electron systems;^ Low- 
energy excitations in a fluid of interacting electrons can 
be treated as excitations of independent quasiparticles, 
whose energy-momentum relationship (the quasiparticle 
energy band) generally differs from that of free electrons. 
Close to the Fermi surface, the quasiparticle band can 
be approximated by the free-particle form appropriate 
for particles of mass m*, where the quasiparticle effec- 
tive mass m* may differ from the bare mass of an elec- 
tron. Given the widespread use of Fermi liquid theory, 
it comes as something of a surprise to learn that the ef- 
fective mass of a paramagnetic two-dimensional (2D) ho- 
mogeneous electron gas (HEG) has been the subject of 
great controversy in recent years. Early experiments^ 
found a large enhancement of the effective mass at low 
density, but subsequent experiments^!^ have found the 
increase in the effective mass to be considerably smaller. 
On the theoretical side, many-body perturbation theory 
{GW) calculations give a range of possible results de- 
pending on the choice of effective interaction and whether 
or not the Dyson equation is solved self-consistently^^i^ 
while quantum Monte Carlo (QMC) studies have found 
either much less^iii^ or much greater— enhancement of 
the effective mass than the experiments suggest. Finally, 
some recent experiments^ have shown that the effective 
masses in paramagnetic and ferromagnetic HEGs behave 
quite differently as a function of density, as had been 
predicted using many-body perturbation theoryJ^ The 
experiments show that the effective mass of a ferromag- 
netic HEG decreases as the density is lowered, which has 
also been observed in a recent GW studyJ^ Understand- 
ing the magnetic behavior of the 2D HEG at low density 
will play an important role in the design of spintronic 
devices. 

In this article, we present QMC calculationsii of the 
energy band of the 2D HEG. We have calculated the 



band £{k) by evaluating the difference in the total en- 
ergy when an electron is added to or removed from a 
particular momentum state k. To our knowledge, this is 
the first QMC calculation of the complete 2D HEG occu- 
pied energy band. As explained above, electronic excita- 
tions close to the Fermi surface correspond to quasipar- 
ticle excitations. The electronic and quasiparticle bands 
therefore agree near the Fermi surface and have the same 
derivative at kp. The effective mass of a HEG can be 
written as^ m* = kp/{d£ /dk)kp, where kp is the Fermi 
wave vector and hence it is straightforward to compute 
the effective mass once the energy band has been deter- 
mined. Our effective-mass data should help to resolve 
the controversies surrounding the increase in the effec- 
tive mass of the paramagnetic 2D HEG at low density. 
We have studied both paramagnetic and ferromagnetic 
HEGs in order to look for the differences in behavior ob- 
served by Padmanabhan et al^ 

We use Hartree atomic units {h = \e\ = rrie = 47reo = 
1) throughout. Densities are given in terms of the radius 
rs of the circle that contains one electron on average. All 
our QMC calculations were performed using the CASINO 
codeJ^ 

The rest of this article is arranged as follows. We ex- 
plain our methodology in Sec. [Hi We present our results 
in Sec. IHIl Finally we draw our conclusions in Sec. |lVl 



II. METHODOLOGY 

In the variational quantum Monte Carlo (VMC) 
method, expectation values are calculated with respect to 
a trial wave function, the integrals being performed by a 
Monte Carlo technique. In the diffusion quantum Monte 
Carloi^ (DMC) method the imaginary-time Schrodinger 
equation is used to evolve an ensemble of electronic con- 
figurations towards the ground state. Fermionic sym- 
metry is maintained by the fixed-node approximation,— 
in which the nodal surface of the wave function is con- 
strained to equal that of a trial wave function. 
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Our trial wave functions consisted of Slater determi- 
nants of plane-wave orbitals multiplied by a Jastrow 
correlation factor. The Jastrow factor contained poly- 
nomial and plane-wave expansions in electron-electron 
separations^ The orbitals in the Slater wave function 
were evaluated at quasiparticle coordinates related to 
the actual electron positions by backflow functions con- 
sisting of polynomial expansions in electron-electron 
separations^ The wave functions were optimized by vari- 
ance minimizatio n^°i ^ ^ followed by linear-least-squares 
energy minimizatioui^ The high quality of our trial wave 
functions is demonstrated in this paper and in Ref. [2^. 

The single-particle energy for an occupied state at 
wave vector k is defined to be £(k) = Eq — £'_(k), 
while the single-particle energy for an unoccupied state 
is f (k) = i5+(k) — Eq, where Eq is the ground-state 
total energy, £'+(k) is the total energy of the system 
with an extra electron placed in orbital exp(ik • r)^ 
and E- (k) is the total energy with an electron removed 
from orbital exp(ik • r). In a finite simulation cell sub- 
ject to periodic boundary conditions, the available states 
{k} fall on the grid of reciprocal-lattice points offset by 
the simulation-cell Bloch vector k^.— The simulation 
cell was left unchanged when electrons were added or re- 
moved. We have confirmed that finite-size biases are neg- 
ligible by carrying out simulations in different cell sizes: 
see Figs. [2] and [3] We have also verified that the band 
values obtained with different simulation-cell Bloch vec- 
tors lie on the same curve as the band values obtained 
with ks = 0. Having determined the energy band at a 
series of k values, we performed a least-squares fit of a 
quartic function £{k) = ao -I- a2k^ + a^k^ to the band 
values. The DMC energy band is defined as a difference 
in total-energy eigenstates; as explained in the introduc- 
tion, this coincides with the quasiparticle band near the 
Fermi surface and hence gives a correct description of the 
effective mass. 

The number of electrons N in each of our ground-state 
calculations was chosen to be a "magic number" corre- 
sponding to a closed-shell configuration when kg = 0. 
In this case real, single-determinant wave functions are 
appropriate for the ground-state calculations, facilitating 
the optimization of the wave function. In the {N + 1)- 
and {N — l)-electron excited-state calculations, we used 
the Jastrow factor and backfiow function that were op- 
timized for the iV-electron ground state. Reoptimizing 
the wave function in the excited state was not found to 
make a significant difference to the VMC or DMC ener- 
gies. To try to obtain a better estimate of the energy 
of the (TV -I- l)-electron system, we constructed a mul- 
tideterminant wave function in which the extra electron 
occupied each of the symmetry-equivalent k vectors in 
the partially filled shell. The determinant coefficients 
were free parameters, which we optimized by linear-least- 
squares energy minimization. However, we were unable 
to lower the VMC energy significantly using this form of 
wave function, so in our production calculations we used 
single-determinant wave functions for the {N + 1)- and 



(iV—1) -electron systems. Our DMC results are converged 
with respect to time step, as is clear from the agreement 
between the energy bands obtained at different time steps 
in Fig.lUb). 

Unlike the DMC calculations of Kwon et al.r-'^ we did 
not promote electrons from the ground-state configura- 
tion; we simply added or subtracted single electrons. The 
energy difference that results from promoting an electron 
contains a contribution from the interaction between the 
excited electron and the hole that it leaves behind, in 
addition to the difference of band energies. By contrast, 
the energy difference resulting from adding or subtract- 
ing an electron simply gives the corresponding band en- 
ergy. We believe our approach to be a simpler procedure 
for calculating the band and hence the effective mass, 
although it does not give values for the quasiparticle- 
interaction Fermi-liquid parameters. Our determination 
of the effective mass will facilitate subsequent calculation 
of the other Fermi-liquid parameters using the approach 
of Kwon et al. 

The other important difference between our methodol- 
ogy and that of Kwon et al. is that we have used a fit to 
the entire occupied band to determine the derivative at 
kp and hence the effective mass. This was done for the 
following reasons: (i) evaluating the derivative numeri- 
cally using only a few band values near fc^? is unreliable 
because of the noise in the band data; (ii) the DMC- 
calculated band suffers from Hartree-Fock-like patholog- 
ical behavior in the vicinity of kp because the method 
does not retrieve all the correlation energy (see the dis- 
cussion in Sec. IIIIBP : and (iii) the band evaluated in a 
finite cell may suffer from finite-size effects in the vicinity 
of the Fermi surfaceJ^ Although the pathological behav- 
ior dominates the derivative of the band in the vicinity of 
kp, it has only a negligible effect on the band fitted over 
a wide range of k. The fit to the band is good, so that 
the derivative of the fitted band at fc^? should be reliable: 
see Figs. [2] and [3] and Figs. [8] and M 

The occupied bandwidth of the HEG is AS ~ £{kp) — 
£{Q) = E^{Q) - E^{kp). The DMC bandwidth is ex- 
pected to be an upper bound on the true bandwidth: 
assuming that DMC retrieves the same fraction of the 
correlation energy in the ground and excited states, the 
bandwidth will lie between the Hartree-Fock (HF) value 
E^^{t)) ~ E^^{kp), which is too large^ and the ex- 
act resuh E°^^'^\0) - E^_^'"^\kp). Likewise, the Slater- 
Jastrow DMC bandwidths are expected to be greater 
than the Slater- Jastrow-backflow DMC bandwidths, as 
can be seen to be the case in Fig. [H^b). To obtain an ac- 
curate bandwidth, it is essential to retrieve a very large 
fraction of the correlation energy in the QMC calcula- 
tions, which explains why the use of DMC and the inclu- 
sion of backfiow is so important in this work. 

A crude way of estimating the ground-state energy is 
to plot the VMC energy against the variance obtained 
with different trial wave functions and extrapolate the 
VMC energy linearly to zero variance, as shown in Fig. 
[T] This procedure suggests that our DMC calculations re- 
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TABLE I: Band widths of paramagnetic and ferromagnetic 2D 
HEGs of density parameter r^, as calculated using DMC, free- 
electron theory {AS = k%/2), and HF theory [AS = k%/2 + 
kpil — 2/7r)]. The DMC bandwidths were obtained from the 
fitted curves shown in Figs. [2] and [3] 



FIG. 1: (Color online) VMC variance against energy for differ- 
ent trial wave functions for a 58-electron paramagnetic HEG 
of density parameter = 5 a.u. Plots are shown for the 
ground state (GS) and an excited state in which an elec- 
tron is removed from the highest-occupied shell. The slanted 
lines show fits to the VMC data. The vertical lines show 
the fixed-node DMC energies obtained with Slater- Jastrow- 
backflow wave functions. The HF ground-state and excited- 
state energies are —0.100222 and —0.099632 a.u. per electron, 
respectively. The difference between the slanted and vertical 
lines at zero variance gives an approximation to the correla- 
tion energy missing in the DMC calculation. This is clearly 
small compared with the difference between the HF and DMC 
energies. 



trieve more than 99% of the correlation energy, and that 
the fraction retrieved is similar in both the ground-state 
and excited-state calculations. Suppose the free-electron 
bandwidth is greater than or approximately equal to the 
exact bandwidth (this is true for the ferromagnetic HEG 
and approximately true for the paramagnetic HEG), so 
that the error in the HF bandwidth is less than or ap- 
proximately equal to A£HF _ ^ghec ^ fc^(i - 2/7r). 
Hence the error in the DMC bandwidth is less than 
0.01fcF(l -2/7r) w 0.007/rs for a ferromagnetic HEG and 
less than about 0.01kp{l — 2/tt) w 0.005/rs for a param- 
agnetic HEG. Since the bandwidth falls off as rj^, the 
error is more significant at large r,,. In the worst case, the 
paramagnetic HEG at rg = 10 a.u., this argument sug- 
gests that DMC overestimates the bandwidth by about 
9%. The errors in the other results are much smaller. 
It seems reasonable to assume that if the bandwidth is 
overestimated by a given amount then the effective mass 
will be underestimated by a similar fraction. 



III. RESULTS 
A. Energy bands 

Our calculated energy bands are shown in Figs. [3] and [3] 
for paramagnetic and ferromagnetic HEGs, respectively. 



The free-electron and HF bands are shown for compar- 
ison. As is well-known, the free-electron band is very 
much more accurate than the HF band, especially at 
low densities and especially in the paramagnetic fluid. 
The HF band is pathological due to the long range of 
the exchange hole, which results in incomplete screen- 
ing of the Coulomb interaction!^ Our DMC bandwidths 
are shown in Table [J The bandwidth in paramagnetic 
HEGs at intermediate and low densities is less than the 
free-electron bandwidth. On the other hand, in the fer- 
romagnetic HEG the bandwidth is greater than the free- 
electron bandwidth at all densities studied. In each case 
the DMC bandwidth is considerably smaller than the HF 
bandwidth. 

It is striking how closely the DMC band agrees with the 
free-electron band £{k) = for a paramagnetic HEG 
at Ts = 1 a.u. At = 5 and 10 a.u., the quartic term 
aik'^ in the fitted band is positive for the paramagnetic 
HEG. For the ferromagnetic HEG the quartic term is 
negative at all densities. In either case the band is clearly 
not quadratic. This will result in non-free-particle-like 
thermodynamic behavior at high temperatures. 



B. Effective masses 

Our DMC effective masses for paramagnetic and ferro- 
magnetic HEGs are plotted in Figs. |4] and [5l respectively, 
along with various experimental results and previous the- 
oretical predictions. Our effective masses are also given 
in Table [Hi In a paramagnetic HEG the effective mass 
increases with rg: at = 1 a.u. the effective mass is 
slightly less than the bare electron mass, but at Ts = 5 
a.u. the effective mass is significantly enhanced. On the 
other hand, in ferromagnetic HEGs the effective mass de- 
creases when the density is lowered. Our results therefore 
support the conclusions of Padmanabhan et al^ In fact 
our ferromagnetic effective masses are in good quantita- 
tive agreement with the experimental data of Padmanab- 
han et al.^ while our paramagnetic effective masses are 
in reasonable agreement with those measured by Tan et 
al^ We do not find especially good agreement with earlier 
theoretical work, however. As can be seen in Fig. HI the 
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FIG. 2: (Color online) Energy bands of paramagnetic A'^- 
electron 2D HEGs at = 1 a.u. (top), 5 a.u. (middle), and 
10 a.u. (bottom). The free-electron and HE bands are off- 
set to coincide with the fitted DMC band at k = kp- The 
curve labeled "SJ" used a Slater- Jastrow trial wave function; 
the others used a Slater-Jastrow-backflow trial wave function. 
Except where indicated otherwise, DMC time steps r of 0.01, 
0.2, and 0.4 a.u. were used at = 1, 5, and 10 a.u. The solid 
lines show quartic fits to the DMC data for N = 74, 114, and 
114 electrons at — 1, 5, and 10 a.u., respectively. (These 
are the largest system sizes for which we have sufficient data 
to perform an adequate fit.) 



EIG. 3: (Color onhne) As Eig.[2l but for ferromagnetic HEGs. 
The solid lines show quartic fits to the DMC data for = 57, 
101, and 57 electrons at Vs — 1, 5, and 10 a.u., respectively. 



effective masses obtained using the GW method depend 
strongly on the choice of effective interaction and whether 
or not the Dyson equation is solved self-consistently, un- 
dermining confidence in that approach. The DMC data 
of Kwon et al. do not show a significant enhancement of 
the paramagnetic effective mass at low densitiesi^ (Our 
calculations differ from those of Kwon et al. in that we 
use wave functions that retrieve a greater fraction of the 
correlation energy, we use larger system sizes, and we use 
different excitations to evaluate the effective inass.) 

Holzmann et alr^ have recently studied the paramag- 
netic 2D HEG effective mass using the VMC method. 
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FIG. 4: (Color online) Effective mass m* against density pa- 
rameter rs for paramagnetic or partially spin-polarized 2D 
HEGs, as calculated or measured by different authors. Our 
DMC results were obtained from the fitted curves shown in 
Fig. [2] The GW results were obtained using the random- 
phase-approximation (RPA) effective interaction^ and the 
Kukkonen-Overhauser (KO) effective interaction^ by solving 
the Dyson equation self-consistently (SC) or within the on- 
shell approximation (OSA). All the results shown are for 
paramagnetic HEGs with the exception of the experimental 
results of Ref.@, which are for a partially spin-polarized HEG. 
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FIG. 5: (Color online) Effective mass m* against density pa- 
rameter rs for ferromagnetic 2D HEGs. Our DMC results 
were obtained from the fitted curves shown in Fig. [S] The 
GW results were obtained using the Kukkonen-Overhauser 
(KO) effective interaction by solving the Dyson equation 
self-consistently (SG) or within the on-shell approximation 
(OSA).^ 
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TABLE II: DMC-calculated effective masses m* of paramag- 
netic and ferromagnetic 2D HEGs of density parameter ra. 



Their cfFcctive masses differ significantly from our DMC 
results: see Fig. S) Holzmann et al. considered additions 
of electrons in the vicinity of the Fermi surface. They 
found substantial finite-size effects in the effective mass, 
which they corrected by considering the finite-size depen- 
dence of the discontinuity in the momentum distribution 



FIG. 6: (Color online) DMC effective mass m* against system 
size A*' for paramagnetic 2D HEGs. 



at the Fermi edge. Our effective masses do not appear 
to suffer from these finite-size effects; in fact, our effec- 
tive masses show the opposite trend with system size, as 
can be seen in Fig. [6l If the finite-size correction to the 
mass falls off as N~'^/^ as predicted by Holzmann et al. 
then the correction should be roughly halved on going 
from iV = 18 to = 202 electrons. Noting that the 
correction is supposed to increase the mass, this ought 
to be visible in our data. Since it is not, the finite-size 
correction of Holzmann et al. appears to be statistically 
insignificant. 

A possible explanation for this difference is that we 
obtained our effective masses by fitting a band to a wide 
range of k values instead of just considering the bchav- 
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FIG. 7: (Color online) HF energy bands (a) and their deriva- 
tives (b) for A^-electron paramagnetic 2D HEGs of density 
parameter = 10 a.u. 



lor in the vicinity of the Fermi surface. As can be seen 
in Fig. [3 the HF bands of finite systems as well as in- 
finite ones exhibit pathological behavior near the Fermi 
surface. The derivatives of the bands become large and 
fluctuate wildly. The numerical derivatives of the DMC 
bands arc plotted in Figs. [5] and [HI It can be seen that 
the DMC bands exhibit residual Hartrec-Fock-likc patho- 
logical behavior at the Fermi surface, because not all the 
correlation energy is retrieved. The pathological behav- 
ior is more pronounced in paramagnetic HEGs and at low 
densities. Our procedure of fltting the band over a wide 
range of k enables us to avoid this pathological behav- 
ior. The finite-size effects considered by Holzmann et al. 
only affect the energy band in the vicinity of the Fermi 
surface, suggesting that it is also possible to avoid these 
finite-size effects by fitting the band over a wide range of 
k. 

It may seem counterintuitive to use DMC data for ex- 
citations far from the Fermi surface in order to determine 
the effective mass, which is a parameter in a theory de- 
scribing excitations in the vicinity of the Fermi surface. 
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FIG. 8: (Color online) Derivatives of the energy bands of 
paramagnetic A''-electron 2D HEGs at = 1 a.u. (top), 5 
a.u. (middle), and 10 a.u. (bottom). The curve labeled "SJ" 
used a Slater- Jastrow trial wave function; the others used a 
Slater- Jastrow-backfiow trial wave function. Except where 
indicated otherwise, DMC time steps r of 0.01, 0.2, and 0.4 
a.u. were used at Ts = 1, 5, and 10 a.u. The central difference 
approximation was used to evaluate the numerical derivatives. 
The solid lines show the derivatives of quartic fits to the DMC 
energy bands for A'^ = 74, 114, and 114 electrons at rs = 1, 5, 
and 10 a.u., respectively. 



However, we reiterate that it is a premise of Fermi liq- 
uid theory that the quasiparticle band coincides with the 
energy band defined by differences in total-energy eigen- 
values in the vicinity of the Fermi surface, so that the 
derivative of the quasiparticle band at the Fermi surface 
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FIG. 9: (Color online) As Fig. [HI but for ferromagnetic HEGs. 
The solid lines show the derivatives of quartic fits to the DMC 
energy bands for N = 57, 101, and 57 electrons at = 1, 5, 
and 10 a.u., respectively. 



is equal to the derivative of the electronic energy band. 
It is possible in principle that the gradient of the energy 
band may change sharply in the vicinity of the Fermi sur- 
face, but there seems to be no reason to suppose this to 
be the case and our results do not provide any evidence 
for this sort of behavior: the energy bands shown in Figs. 
[5] and [3] look well-behaved. Nevertheless, it is comforting 
to observe that our effective-mass results are insensitive 
to a reduction in the range of fc-values used to perform 
the fit. This result is obvious from looking at Figs. [2] 
and [21 and it also follows from a quantitative study, as 
shown in Tables IIIII and IIVI For both the paramagnetic 



Range 




No. pts in fit 


m* (a.u.) 


0.0 < k/kp 


< 1.33 


17 


1.21(2) 
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< 1.26 
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1.18(2) 
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< 1.20 


13 


1.14(2) 


0.47 < k/kp 


< 1.17 


11 


1.12(3) 



TABLE III: Effective mass m* versus range of k values used 
to fit the energy band for a 114-electron paramagnetic HEG 
at r-s = 5 a.u. 



Range 


No. pts in fit 


m* (a.u.) 


0.0 < k/kp < 1.18 


12 


0.706(9) 


0.18 < k/kp < 1.12 


10 


0.72(1) 


0.35 < k/kp < 1.06 


8 


0.722(9) 


0.5 < k/kp < 0.95 


6 


0.686(8) 



TABLE IV: Effective mass m* versus range of k values used 
to fit the energy band for a 101-electron ferromagnetic HEG 
ai rs ~ 5 a.u. 



and ferromagnetic HEGs, discarding data far from kp 
has little effect on the calculated effective mass, until the 
remaining data points are all sufficiently close to kp that 
either the pathological behavior of the DMC energy band 
in this region starts to dominate or there are insufficient 
data to perform an accurate fit. 

Another possible reason for not reproducing the finite- 
size errors in the effective mass predicted by Holzmann et 
al. might simply be that some of their assumptions are in- 
valid. We have reproduced their 0{N~^^^) scaling of the 
finite-size error in the renormalization factor (the discon- 
tinuity Z in the momentum density at the Fermi edge) in 
VMC calculations, as shown in Fig. [TOl However, when 
an electron is added, there arc two contributions to the 
momentum density: a peak of weight Z at the momen- 
tum at which the electron is added and a smeared-out 
background of weight 1 — Z.— Together, these two con- 
tributions to the change in the momentum density are 
responsible for the change in the kinetic energy when 
the electron is added, i.e., for the kinetic contribution to 
the energy band. Examples of the changes in the VMC 
momentum density that result from adding or removing 
electrons from different k are shown in Fig.[TTJ It can be 
seen that the smeared-out background depends on the k 
at which the electron is added, its average tending to in- 
crease with k. The finite-size error in the weight Z of the 
peak is equal and opposite to the finite-size error in the 
weight 1 — Z of the background. The finite-size error in 
the derivative of the energy band due to the background 
therfore tends to cancel the finite-size error due to the 
peak. So it is not clear that the 0{N~^/'^) finite-size 
error in the renormalization factor should result in an 
0{N~^^^) error in the gradient of the energy band and 
hence effective mass. 
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FIG. 10: (Color online) VMC renormalization factor Z 
against number of electrons A'' for paramagnetic 2D HEGs 
at Ts — 10 a.u. 



0.04 



0.02 



-0.02 




(Addition, 1st shell 
I Addition, 2nd shell 
► Addition, 3rd shell 
k Subtraction, 1st shell 
Subtraction, 2nd shell 
Subtraction, k = 
I 



k/k^ 



2.5 



FIG. 11: (Color online) VMC momentum density p{k) relative 
to the ground-state momentum density pGs(fc) for different 
excitations to a paramagnetic 58-electron HEG at Vs = 10 
a.u. The momentum densities are averaged over reciprocal- 
lattice vectors of the same length, so the height of the spike 
at the |k| at which the electron is added looks smaller when 
there are many reciprocal-lattice vectors with the same length 
as Ikl. 



In summary, we have used DMC to calculate the en- 
ergy band of the interacting 2D HEG, and hence we have 
obtained the quasiparticle effective mass. Our ferromag- 
netic and paramagnetic effective masses are in agreement 
with the experimental results of Padmanabhan et al^ 
and Tan et al.^ respectively. In particular, our data 
confirm that the effective mass of the paramagnetic HEG 
increases when the density is lowered, while the effective 
mass of the ferromagnetic HEG decreases. 
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